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The entanglement properties of a class of topological stabilizer states, the so called topological 
color codes defined on a two-dimensional lattice or 2-colex, are calculated. The topological entropy 
is used to measure the entanglement of different bipartitions of the 2-colex. The dependency of 
the ground state degeneracy on the genus of the surface shows that the color code can support a 
topological order, and the contribution of the color in its structure makes it interesting to compare 
with the Kitaev's toric code. While a qubit is maximally entangled with rest of the system, two 
qubits are no longer entangled showing that the color code is genuinely multipartite entangled. For 
a convex region, it is found that entanglement entropy depends only on the degrees of freedom 
living on the boundary of two subsystems. The boundary scaling of entropy is supplemented with 
a topological subleading term which for a color code defined on a compact surface is twice than the 
toric code. From the entanglement entropy we construct a set of bipartitions in which the diverging 
term arising from the boundary term is washed out, and the remaining non- vanishing term will 
have a topological nature. Besides the color code on the compact surface, we also analyze the 
entanglement properties of a version of color code with border, i.e triangular color code. 

PACS numbers: 03.67.Lx, 03.67.Mn, 03.65.Ud, 42.50.Dv 
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INTRODUCTION 



In quantum information theory and quantum compu- 
tations, entanglement is recognized as an essential re- 
source for quantum processing and quantum communi- 
cations, and it is believed that the protocols based on 
the entangled states have an exponential speed-up than 
the classical ones. Besides, in highly correlated states 
in condensed matter systems such as superconductorsi*^, 
fractional quantum Hall liquids^, the entanglement serves 
as a unique measure of quantum correlations between de- 
grees of freedom. In quantum many-body systems such 
as spin, fermion and boson systems, the entanglement is 
connected to the phase diagram of the physical systems. 
Non-analytic behavior of entanglement close to the quan- 
tum critical point of the system and occurrence of finite 
size scaling have provided an intense research leading to 
fresh our insight of the critical properties from the quan- 
tum information side, for a comprehensive discussion see 
the review by L. Amico, et alA- and references therein. 
In the past years the appearance of new phases of mat- 
ter has intensified the investigation of the entanglement 
in the quantum systems. These are phases beyond the 
Landau-Ginzburg- Wilson paradigm^ where an appropri- 
ate local order parameter characterizes different behav- 
iors of two phases on either side of the critical point. 
These new phases of matter carries a kind of quantum 
order called topological order^ and the transition among 
various phases does not depend on the symmetry break- 
ing mechanism. Therefore the Landau theory of classical 
phase transition fails in order to describe these phases. 
Ground state of such phases is a highly entangled state 
and the excitations above the ground state have a topo- 
logical nature which mirrored in their exotic statistics. 

Among the models with topological properties, the Ki- 
taev or toric code^ has been extensively studied. Ground 



state degeneracy depends on the genus or handles of the 
manifold where the model is defined on, and there is a 
gap which separated the ground state subspace from the 
excited states. The ground state degeneracy can not be 
lifted by any local perturbations which underlines the 
Kitaev's model as a testground for fault-tolerant quan- 
tum computations^. The ground state of the Kitaev's 
model is indeed stabilized by group generated by a set 
of local operators called plaquette and star operators 
making it useful as a quantum error correcting code^. 
The information is encoded in the ground state subspace 
which is topologically protected making it robust as a 
quantum memory^. For this model any bipartition of 
the lattice has non-zero entanglement which manifests 
the ground state is generically multipartite entanglediS. 
Bipartite entanglement scales with the boundary of the 
subsystem showing that the entanglement between two 
parts of the system depend only on the degrees of free- 
dom living on the boundary which is a manifestation of 
the holographic character of the entanglement entropyii. 

For topological models a satisfactory connection 
between topological order and entanglement content of a 
model has been established via introducing the concept 
of the topological entanglement entropy (TEE)^^^ 
which is a universal quantity with topological nature. 

Another resource with topological protection charac- 
ter is called color code. An interplay between color and 
homology provides some essential features, for example 
a particular class of two-dimensional color codes with 
colored borders will suppress the need for selective ad- 
dressing to qubits through implementation of Clifford 
group^'^.The number of logical qubits which are encoded 
by a two dimensional color code are twice than the toric 
code^^defined on the compact surface. 
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In this paper we study the entanglement properties of a 
color code defined on a two-dimensional lattice, the so 
called 2-colex— . We consider different bipartitions of the 
lattice and evaluate the entanglement entropy between 
them. The connection between the topological nature of 
the code and the entanglement entropy is also discussed. 
The structure of the paper is as follows. In sec (II) the 
basic notions of a color code on the surface has been in- 
troduced. In sec(III) based on the stabilizer structure 
of the protected code space the reduced density matrix 
which is needed for evaluating the entanglement entropy 
has been calculated. Then in sec(IV) the entanglement 
entropy for different bipartitions is calculated. In sec(V) 
another class of color code on plane and its entanglement 
properties are introduced, and the final section (VI) has 
been devoted to the conclusions. 



II. PRELIMINARIES ON TOPOLOGICAL 
COLOR CODE 

In this paper we consider a class of topological quan- 
tum error-correction code defined on the lattice, the so 
called topological color codes (TCC)i^. The local degrees 
of freedom are spin- 1/2 with the bases of the Hilbert 
space C^. The lattice we consider composed of vertices, 
links and plaquettes. Each vertex stands for a local 
spin. Three links meet each other at a vertex and no 
two plaquettes with same color share the same link. We 
suppose this color structure is denoted by the notation 
TCC{V,E,P} where the V,E,P denote the set of ver- 
tices, edges and plaquettes, respectively. For simplicity 
we define the model on the regular hexagonal lattice on 
the torus, i.e imposing periodic boundary conditions as 
shown in Fig.([T]). There is a subspace C C H which is 
topologically protected. The full structure of this sub- 
space and its properties are determined by definition of 
stabilizer group. The stabilizer group is generated by a 
set of plaquette operators. For each plaquette we attach 
the following operators which are product of a set of Pauli 
operators of vertices around a plaquette: 



(^rj^; Q^X,Z, C = Red, Green, Blue. (1) 



For a generic plaquette, say blue plaquette Pi in 
Fig-lHJ we can identify green and red strings which are 
the boundary of the plaquette. It is natural to think of 
product of different plaquette operators which may pro- 
duce a collection of boundary operators. For example as 
is shown in Fig.([T]) the product of two neighboring pla- 
quettes, say red and blue ones, correspond to a green 
string P2 which is a boundary of two plaquettes. All 
string operators produced in this way are closed. Since 
all closed strings share either nothing or even number of 
vertices, they commute with each other and with plaque- 
ttes . In addition to closed boundary operators, there 
are other closed string which are no longer the product 
of the plaquette operators. These closed string form the 



fundamental cycles of the manifold in which the lattice is 
defined on and have a character of color. Number of these 
closed loops depends on the genus of the manifold where 
the lattice defined on. For the torus with g = 1 there 
are two such cycles which are non-contractible loops in 
contrary to the closed boundary strings which are homo- 
topic to the boundary of a plaquette. For the topologi- 
cal color codes these non-contractible loops are shown in 
Fig. (II]). For every homology class of the torus there are 
two closed strings each of one color, say red and blue. 
Red string connects red plaquettes and so on. Note that 
a generic string, say green, can be produced by product 
of the red and blue strings when they are suitably cho- 
sen. In fact since every Pauli matrix squares identity, at 
a vertex (qubit) which two homologous strings cross each 
other they cancel each other. Indeed there is an interplay 
between color and homology class of the model. One can 
define a nontrivial closed string as follows: 



5f ^ = ' 



(2) 



where / indexed the set of spins on a generic string, fi 
stands for the homology class of the torus and fl is the X 
or Z Pauli spin operators. Closed non-contractible loops 
turn on to form bases for the encoded logical operators of 
topological code. To more clarify this, we label different 
loops as: 



Zi 



Si 



BX 



qBX y 
O2 7 ^4 



St 



BZ 



Si 



RZ 



Si 



RX 



.(3) 



These operators form a 4-qubit algebra in 7^1, so it 
manifests a 16-dimensional subspace for the coding space 
C which is topologically protected. On the other hand, 
the topological color code on the torus with g = 1 encodes 
four qubits. Now we move to construct the explicit form 
of the states of the subspace C. 

The above construction for the string operators in Eq. ([3]) 
can be extended to an arbitrary manifold with genus g. 
For such manifold the coding space spanned with 2^^ 
vectors. Note that for toric code (white and dark code)^^ 
when is embedded in the same manifold with genus g, the 
coding space will span with 2^^ vectors^*' which explicitly 
shows that the color codes have richer structure than the 
toric codesi^. 



A. Stabilizer Formalism 

The protected subspace C is spanned by the state 
vectors which are stabilized by all element of stabilizer 
group, i.e. a subset of Pauli group. Let W be a set of 
generators of the stabilizer group and its elements are 
denoted by A^. So this subspace is 
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FIG. 1: (color online) A piece of 2-colex which has been de- 
fined on the hexagonal lattice with periodic boundary con- 
ditions. Lightest to darkest hexagons: green, red and blue, 
respectively. All circles stand for qubits living at the ver- 
tices. A hexagon (Pi), say blue, can be denoted by either 
red or green closed boundary string and the product of two 
neighboring plaquettes, say red and blue, corresponds to a 
green string (^2). Non-contractible loops {Xi,X2,X3,X4,) for 
topological color codes which determine the structure of the 
encoded subspace C. A colored string connect plaquettes with 
the same color. 



C = {IV) : - 10) VMeU} 



(4) 



Let G be the group constructed by the generators of 
spin-flip plaquette operators, i.e X^ . The cardinaUty of 
the group is then \G\ = 2^^^'^ where \P\ stands for the 
total number of plaquettes. Note that all plaquettes are 
not independent since the product of all plaquettes with 
the same color represents the same action in the group, 
namely U^ Xp = Ub = Ug^p = where 
\V\ stands for the number of all vertices. By starting from 
an initial vacuum state, say jO)®'^' where Z\0) — |0), 
one can construct a state vector in the Hilbert space 
which is stabilized by the group elements. This state 
vector is a superposition of all elements of the stabilizer 
group with equal weights. From the previous arguments, 
it is convenient to denote it by |0000) which has the fol- 
lowing form 



10000) = |G| 



-/2^,|0)C 



(5) 



where the g is an element of the stabilizing group, i.e 



5 = 0„eP^P 



where Vp — 0(1) corresponds to the pla- 



quette operator Xp appearing (not appearing) in the el- 
ement group g. There are many elements of the Pauli 
group that commute with the all elements of the stabi- 
lizer group but are not actually in G, and it is defined 



as the centralizer of G in Pauli group. Since elements 
of Pauli group either commute or anticommute, the cen- 
tralizer is actually equal to the normalizer of G in Pauli 
group. Considering the non-contractible loops, the nor- 
malizer of stabilizer group can be obtained by product of 
stabilizer group and group of non-contractible loop op- 
erators. Let denote it by G = .4 • G where the group A 
generated by non-contractible loop operators, i.e. 



-4={ n (^i 



CX\ 

fl ) 



0,1} 



(6) 



fj,=l,C=B,R 



So, the group G reads 

G = { n' (S^^'P-'G , r^c = 0,1} . (7) 

Note that G G. The group G is the normalizer of 
G in Pauli group, so the normal subgroup G divides the 
group G into 2^^^ cosets. The cardinality of group will 
be |G| = 21^1+'^^^^. Therefore the protected subspace 
C is spanned by states which correspond to different 
cosets. In fact, elements in G — G take one encoded 
state of the stabilized subspace to another encoded state 
without leaving the stabilized subspace. For the states 
of the stabilized subspace we have, by construction: 

C^{\ijkl): \ijkl) = XlX^X^Xi\0000) }, (8) 

where Xi, X2, ^3, are defined in Eq.Q and i, j,k,l — 
0, 1. These topological non-trivial string operators can 
take one state of the coding space to another one, and 
any error of this type will not be detectable. A generic 
state can be a superposition of different vectors of coding 
space as follows: 



l*) = 



i,j,k,l 



'i,j,k,l 



\ijkl) 



E 

i,j,k,l 



'i,j,k,l I 



= 1 



(9) 



B. Protected Subspace as a Ground State Subspace 

From the practical point of view it is important to find 
a state of quantum many-body system for implementa- 
tion of universal quantum computation. We can provide 
a construction in which the protected subspace be the 
ground sate of a local Hamiltonian^^. The subspace C is 
the ground state of a following exactly solvable Hamilto- 
nian , i.e. 



H - -^Xp-^Zp 

peP peP 



(10) 



The ground state of this Hamiltonian is 2'^S-io\d degen- 
erate and topologically protected from local errors. Dif- 
ferent states of the ground state subspace as it is clear 
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from Eq. ([8]) are obtained by product of non-contractible 
colored strings. Each state is characterized by a set 
of topological numbers which are sum of all z^, mod- 
ula 2 along non-contractible loops. For example for a 
generic state \iikl) in Eq.® the topological numbers 
are E«e/f = E^e/« = 3-, E^e/f ^ k and 
Sue/f '^'^ = where the summations are evaluated in 
mod 2, Zv stands for the ^-component of Z-Pauli matrix 
of vertex v and identifies a set of vertices (qubits) on 
a closed non-contractible loop with homology /i and color 
C . An excited state will arise when any of the stabilizing 
conditions in Eq. ([4]) is violated. In the energy unit which 
is defined by Hamiltonian Eq. (fTO|) . the first excited state 
will be separated from the ground state by a gap of value 
2. 



III. REDUCED DENSITY MATRIX AND VON 
NEUMANN ENTROPY 

In this section we turn on to calculate the entangle- 
ment properties of the topological color codes. We con- 
sider a generic bipartition of the system into subsystems 
A and B. Let T,a and to be the number of plaque- 
tte operators acting solely on A and B, respectively, and 
let Yjab to stand for the number of plaquette operators 
acting simultaneously on A and _B, i.e these are bound- 



ary operators. We focus on the entanglement entropy 
between two partitions A and B of the system. To this 
end, first the reduced density operator of the one subsys- 
tem is evaluated and then the entanglement entropy is 
measured. If the state of the system is in an equal super- 
position of the elements of group C?, the reduced density 
matrix for a subsystem, say A, has the following formic 



51 9a\Qa){aQ\9a9a (11) 

where Ga and Gb are subgroups of G which act trivially 
on subsystems B and j4, respectively and cIa and are 
their cardinality. So the von Neumann entropy is 



SA = \0g2\GAB\ (12) 

where Gab = q ■ It is a simple task to show that all 
states of the cocfing space have the same entanglement 
entropy To show this, let X{t) = XlX^X^Xl in which 
t = (i,j,k,l) is a binary vector. So, for a generic state 
in the coding space we have: \t) — X{t)\0). Moreover 
the string operators X{t) can be decomposed as X{t) = 
X{t)A<S> X{t)B- Therefore we obtain 



PA{\t)) = TrBi\t){t\) = TrBiX{t)\0){0\X{t)) = X(<)^Trs(|0)(0|)X(t)^ = X{t)APAm)X{t)A . (13) 

I 



This implies that, by using the ancilla SA(t) = 
lim„^i 9„Tr[p^], the entanglement entropy for a state 
\t) is SA{t) = Sa{0). 

IV. ENTANGLEMENT OF TCC FOR VARIOUS 
BIPARTITIONS 

In this section we design different spin configurations 
as subsystems and then evaluate their entanglement with 
its complementary. 

A. One Spin 

As first example we consider the entanglement between 
one spin and remaining ones of the lattice. In this case 
there is no closed boundary operators acting exclusively 
on the spin, i.e subsystem A. So the reduced density ma- 
trix PA is diagonal as follows: 

PA = .r' 9a\Oa){aO\9a (14) 

g&G/GB 



where / — \G/Gb\ is the number of operators which 
act freely on the subsystem A. As it is clear there are 
three plaquette operators which act freely on one spin 
since every spin in the color code is shared by the three 
plaquettes (See Fig.©). This leads to / = 2^ = 8. Only 
half of operators of quotient group leading to flip the spin, 
i.e. they will have trivial effect unless three plaquettes act 
either individually or altogether. Therefore the reduced 
density matrix is 



Pa = ^{\1a){aI\ + \Oa){aO\). (15) 

Thus, in topological color code every spin is maximally 
entangled with other spins. For the surface code (Kitaev 
model) a single spin is also maximally entangled with 
other spinsiS. 

B. Two Spins 

We calculate the entanglement between two spins. For 
this case also there is not any closed boundary operator 
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with nontrivial effect on two spins, so the reduced density 
matrix reads 



PA = ^(|11)(11| + |10)(10| + |01>(01| + |00)(00|). (16) 

The entanglement between two spins which is mea- 
sured by the Concurrence^* vanishes. While each spin 
is maximally entangled with others, two spins are not 
entangled which is a manifestation of the fact that the 
topological color code is genuinely multipartite entangled 
like the surface codes. 



C. Colored Spin Chain 

In this case we aim to know how much a closed colored 
spin chain winding the torus nontrivially, say blue or red 
in Fig. IT]) is entangled with rest of the lattice. For the 
sake of clarity and without loss of generality we consider 
a hexagonal lattice with \P\ = 3k x 4k plaquettes, where 
k is an integer number as 1,2,3..., for example in Fig.([T]) 
k = 2. This choice makes the color code more symmetric. 
So, the number of plaquettes with a specific color, say 
red, is (2fc)^, and the number of spins that the colored 
chain contains is 4/c. Let the system be in the |0000). 
For this state there is not any closed boundary which 
exclusively acts on the chain. Therefor the reduce density 
matrix is diagonal and the number of plaquette operators 
which act independently on the subsystem B is 

I]s = (|P|-2)-(3x2fc) + (2/c + l) (17) 

where the third term is the number of constrains on the 
plaquette operators acting on spin chain. In fact these 
are collective operators, i.e product of boundary plaque- 
tte operators between A and i?which act solely on B. 
With these remarks in hand, the entanglement entropy 
becomes 



5a = log ^ = 4fc - 1 . (18) 

It is instructive to compare the obtained entanglement 
entropy in Ea. p8)l with the entanglement of the spin 
chain in the Kitaev's model. For the latter case when 
the model defined on the square lattice on a compact sur- 
face, i.e torus, there exists a rather similar scaling with 
the number of qubits living in the chain^". 

In both cases the entanglement entropy in Eq. (fT8|) can 
be understood from the fact that only configurations with 
even number of spin flips of the chain allow in the ground 
state structure. Plaquettes which are free to act on the 



spin chain only able to flip even number of spins of the 
chain. The same arguments can be applied for other 
ground states \ijkl) where for all of them the entangle- 
ment entropy is given by Eq. ljlSp . By inspection of other 
ground states we see that there are 2^ = 8 state in which 
even number of spins in the spin chain have been flipped. 
These states are 

10000) ,Xi|0000) ,X3|0000) ,X4|0000) ,XiX3|0000), 

X1X4IOGOO) , X3X4IOOGO) , 10000). (19) 

Besides, there are 2'^ = 8 states with an odd number 
of spin flipped for the spin chain which have been listed 
below: 

X2IOOOO) , X2X1IOOOO) , X2X3IOOOO), 



X2X4IOOOO) ,X2XiX3|0000) ,X2XiX4|0000), 



X2X3X4IOOOO) ,X2XiX3X4|0000). (20) 

Now we consider a generic state and calculate the en- 
tanglement entropy between the colored spin chain and 
rest of the lattice. The density matrix p = is 

P = 0,ijkiamnpqX\X^X^X\^P()X"^ X2 X^Xl. [21) 

ijkl,mn'pq 

With the explanations we made above the eigenvalues 
of the reduced density matrix pA for the even number of 
spin flipped is as follows: 

/I |2,| |2|| |2|| |2i| |2 

7777 U "^^oooo I + Fioool + loooiol + 1^0001 1 + |aioio| 



+ laiooiT + laooiiT + |aioii|^) = 7777a (22) 

l<^l 

where ^ = 2'^^~^. Note that the number of these 

eigenvalues is 2'^^~^ . Through the similar arguments we 
see that the eigenvalues of reduced density matrix of the 
spin chain with odd number of spin flipped is Jq| (1 ~ ck) 
and the number of these eigenvalues is 2'"''"-'^. So, we can 
calculate the entanglement entropy as follows: 
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Sa^-J2>^^ log ^' = -2''"' log ^« - 2 



where H{x) — —xlogx — (1 — a;)log(l — x). 
What about the entanglement of an open string? the 
open strings in the color code stand for the errors and 
map the coding space into an orthogonal one. These open 
strings anticommute with plaquette operators which 
share in odd number of vertices. For an open string 
which contains k' qubits the entanglement entropy is: 
Sa — log ^ J = fc', i.e an open string will be maximally 
entangled with rest of the system. 

D. Red Strings Crossing 

As an another bipartition we select all spins on the two 
red strings with different homologies, say Xi and X4 in 
Fig. (Ill), which have 2 x 4fc = 8fc spins. Let system be 
in the |0000) state. There is not any closed boundary 
operator acting solely on A. Considering the total num- 
ber of independent plaquette operators acting solely on 
subsystem 5, the entanglement entropy then reads 

SA^Sk^ 1. (24) 

Again we see that there are only configurations with 
an even number of spin flipped in the construction of the 
state 1 0000). For a generic state as in Eq. (PT|) where we 
can realized all states with either even or odd number of 
spin flips, the entanglement entropy will be 

S'a = 8fc - 1 + H{a) . (25) 

E. Red and Blue strings crossing 

By this partition we mean the spins on the Xi and 
X2 strings in Fig.(IT]). The subsystem A contains 8k — 1 
spins and again there is no closed boundary acting on it. 
Enumerating the independent plaquette operators acting 
on the subsystem B, the entanglement entropy is 

S'a = 8fc - 3 . (26) 

F. Two Parallel Spin Chains 

Two parallel spin chains, say red and blue (Xi and X3) 
in Fig.(IT]) or equivalently a green string, contain 2x 4k = 
8k spins and again there is not any closed string acting on 




FIG. 2: (color online) A set of vertical spin ladders (black solid 
circles) on 2-colex as a subsystem for calculating entanglement 
entropy. 



it nontrivially. 8fc plaquettes act freely on subsystem A, 
but there are some constraints on them which arise from 
the product of plaquettes in which leave the subsystem A 
invariant and makes collective operator which acts solely 
on B. The product of blue and green plaquettes of those 
plaquettes which are free to act on A will produce a red 
string leaving A invariant. The same product hold for the 
red and green plaquettes. So the entanglement entropy 
reads 

SA^8k~2 . (27) 

Indeed this is entanglement entropy for a green string. 
However, the number of spins it contains is twice than 
the strings we discussed in Subsec(C). 

G. Spin Ladder 

A set of vertical spin ladders has been shown in Fig. ([2]). 
In this subsection only one of them has been considered 
as subsystem A. Let the system be in the state |0000). In 
this case the subsystem A contains 2 x 3k spins, and the 
total number on independent plaquette acting on sub- 
system i? is (|P| — 2) — 3 X 3fc -|- 3/c. Since there are not 
any closed string acting solely on A, the entanglement 
entropy reads 

SA = Gk . (28) 
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FIG. 3: (color online) A manifestation of a hexagonal disk as 
a subsystem. A set of plaquettes is chosen in such a way that 
form a large hexagon (the yellow line), and all large yellow 
(light) spins are in the subsystem . The dashed black line 
depicts the set of plaquettes acting simultaneously on both 
subsystems. 



So, it is clear that this spin ladder has maximum en- 
tanglement with the rest of the system and its state has 
the maximum mixing, namely 

PA = 2 ®'^l6fcx6fe ■ (29) 

H. Vertical Spin Ladders 

Now we consider all vertical spin ladders has been de- 
picted in Fig.(l2|) and again we suppose the system is in 
the state |0000). There is not any plaquette operator 
which acts solely on subsystem A. However, there are 
several specific product of plaquette operators producing 
closed strings in which act solely on A. For example the 
product of plaquettes which are in a vertical column will 
produce a closed string acting on one of the subsystems. 
Both subsystems A and B are symmetric with respect to 
each other. So, the number of closed strings which act 
only on A and B will he (Ia — ds = 2^*-'. Finally the 
entanglement entropy reads 

Sa = 12P -Ak-2 . (30) 



I. A Hexagonal Disk 

In this case we adapt a situation in which a set of pla- 
quettes intuitively forms a hexagon as shown in Fig. ([3]). 
In this figure the dashed hexagon demonstrates a set 
of plaquettes lying between two subsystems. We sup- 
pose that the number of plaquettes crossed by the one 



edge of the hexagon be n. Now we can enumerate the 
number of plaquettes which act on A in terms of n, i.e. 
YiA = 3n(n — 1) + 1. The total number of spins of the 
subsystem A is 6?!^. There are Y^ab ~ Sn plaquettes 
which act between A and B. Let the system be in the 
state 1 0000). With these realizations of the subsystems 
and since \P\ — "Ea + + ^ab we can turn on to calcu- 
late the entanglement entropy between two subsystems 
as follows: 

S'a = log22^-*«~2 ^ Eab - 2 = 6n- 2 . (31) 

We can relate this entropy to the perimeter of the sub- 
system A by choosing the hight of a plaquette as unit. 
The perimeter of A is dA = 6n and the entropy then will 
be 

SA^dA-2 . (32) 

Generally for a irregular lattice the entropy for a convex 
shape will be Sa = kOA — 7 where the coefficient k is 
a nonuniversal constant depending on the shape of the 
region while the constant 7 = 2 will be universal and 
has a topological nature. So, we see that the entangle- 
ment entropy for the topological color code scales with 
the boundary of the subsystem which is a manifestation 
of so called area la w^^^^^^'^^ . 

The latter relation for we obtained for the entanglement 
entropy is consistent with derivation of A. Kitaev et aZ.— 
and M. Levin et alJ^ where they proposed for a mas- 
sive topological phase there is a topological subleading 
term for the entanglement entropy which is related to the 
quantum dimension of the abelian quasiparticles. The to- 
tal quantum dimension for the topological color code will 
be 

7 = log 2? , D = 4 . (33) 

The quantity is the number of topological superse- 
lection sectors of abelian anyons. For Kitaev's toric code 
there are only four such sectors leading to 7 = log 2. 
So the total quantum dimension in the topological color 
code is bigger than the toric code. The abelian phase 
of the toric code is characterized via appearing a global 
phase for the wavefunction of the system by winding an 
electric (magnetic) excitation around a magnetic (elec- 
tric) excitatiorti^. In the case of color code the excita- 
tions are colored and they appear as end points of open 
colored strings of a shrunk latticelS. Winding a colored 
AT-type excitation around the Z-type one with different 
color gives an overall factor (-) for the wavefunction of 
the model, i.e the phase is abelian. Therefore the contri- 
bution of the color to the excitations makes the abelian 
phase of the color code richer than the toric code with a 
bigger quantum dimension. 

Now let the system be in a generic state such as 
Eg. ipT]) . The reduced density matrix then reads 
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Pa — 2^ O'ijklO.mnpqTrB 
ijkl^mnpq 



Some remarks are in order. For two different nontrivial 
closed string operators but with the same homology and 
color , say X and X' , they will have same support and 
are related via a trivial closed string g E G, i.e X' — gX. 
Since [g,X] — 0, two string will have the same effect on 
the po, i.e. X'po — gXpQ = Xpo- We can exploit this 
property of string operators in order to choose the strings 
appearing in Eq. (|34p in which they do not cross subsys- 



TrB{Xipo)= 9A\0A){AQ\g'AY.(MX4gB\0B){ 
g,9'eG g'^ 




FIG. 4: (color online) Different bipartitions of the lattice in 
order to drop the bulk and boundary effect through the defi- 
nition of topological entanglement entropy. 



On the other hand from the group property g' = gg 
the above expression becomes 



TreiXipo)^ ^ gA\OA) {AO\gAgA{BO\gBXB\OB) (36) 

which implies that gB = Xb- Otherwise the above 
expression will become zero. The fact g = gA ® gB £ G 
explicitly implies that the operator gA must be a non- 
contractile strings which is impossible since it must act 
solely on A. So we will leave with TrsiXipo) = 0. The 
same arguments may also apply for other cases. Finally 
the entanglement for a generic state of the system and a 
convex region will become 



PA=TrB{po) , SA^dA-2 . (37) 



IXiX^XipoX^X^X^X!) . (34) 

I 

tem A. Indeed the strings Xi act only on the subsystem 
B, i.e X,A = Ia- This leads to TreiX.poXi) = TrB{po) 
for i = 1,2,3,4 and TrB{XiXlX^X{poXlXiX^X\) = 
TrB{po) for i,j,k,l — 0,1. Other cases will be zero. 
To clarify this, let for example consider the simple case 
TrB{Xipo) where we have 



0\9'b9b\0b)= 9A\0A){AQ\g'A{B0\g'BXBgB\0B)- (35) 
I 




FIG. 5: (color online) One of subsystems, say A, composed of 
two disjoint regions which introduce some nontrivial colored 
closed strings acting solely on B. The closed blue and red 
strings have been shown only for one region of subsystem A. 



J. Topological Entanglement Entropy 

As we see from the Eg. ([5^. in the entanglement en- 
tropy of a region there exists a subleading term which is 
universal and independent of the shape of region. This is 
a characteristic signature of topological order which has 
been found on the subleading term of the entanglement 
entropy. To be more precise on the topological character 
of subleading term, we can construct a set of bipartitions 



in which the boundary contribution is dropped and re- 
maining term is stemmed from the topological nature of 
the code and inspires the topological order. To drop the 
bulk and boundary degrees of freedom we consider a set 
of partitions which have been shown in Fig.Q. Topolog- 
ical entanglement entropy then arises from the following 
combination of entanglement entropy 

Stopo = lim {-SiA + S2A + S3A ~ S4a) ■ (38) 

R.r — *oo 

For each bipartition one can use the entanglement en- 
tropy as before, namely Sq = log(|G|) — log(c?yid_B)- So 
far the cardinalities cIa and ds were determined by the 
number of plaquette operators acting solely on A and B, 
respectively, i.e. dA = 2^^ and ds = 2^^, where and 
stand for the number of plaquette operators acting 
only on A and B, respectively. However, for calculat- 
ing of the Stopo we should take into account other closed 
strings in addition to the above closed strings acting on 
subsystems. For the case in which each partition is a sin- 
gle connected region, for example bipartitions (2) and (3) 
in Fig. (Ill), the previous arguments work, but for disjoint 
regions such as (1) and (4) in Fig.(|3]) we should extend 
the above result. As an example a simple case has been 
shown in Fig. ([5]) where we have supposed the subsystem 
A is composed of two disjoint regions (two dashed rect- 
angles which we label them as Ai and A2 ) while the 
subsystem i? is a single connected one. By inspection we 
see that the product of two set of plaquettes, say blue 
and green ones, which act on ^1 and the blue and green 
ones acting simultaneously on Ai and B result in a red 
string which acts only on B. The same scenario can be 
applied for other choices of the plaquettes, e.g. red and 
green or red and blue plaquettes, which yield blue and 
green string leaving the subsystem Ai. However, as we 
pointed out before one of them will be immaterial since 
there is an interplay between homology and color. 
The main point is that it is not possible to produce these 
colored closed strings which act only on B from the prod- 
uct of some plaquettes of region B. The same thing will 
be held for the another disjoint region of A, say ^2- So, 
we have a collection of closed strings with nontrivial effect 
on B and we must take into account them in calculation 
oi ds- But, a remark is in order and that if for exam- 
ple we combine two red strings, the resultant is not a 
new closed string because we can produce these two red 
string from the product of the blue and green plaquettes 
of region B. Therefore for the case shown in Fig. ([5]) there 
are only two independent closed strings of this type with 
a nontrivial effect on B. With these remarks, now the 
cardinalities dA and ds are 

dA = 2^-* , dB = 2^«+''-2 . (39) 

We can generalize the above results for the case that 
each subsystem is composed of several disconnected re- 
gions. Let the partitions A and B are composed of ruA 
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(a) (b) 

FIG. 6: (color online) A class of color code on the plane. It 
includes borders each of one color. It encodes one qubits and 
can be produced from a 2-colex without border by removing 
a couple of plaquettes. To each border only strings with the 
same color of border can have end points, (a) 3-string net 
commutes with all plaquette operators. Such string-net and 
its deformation makes the encoded Pauli operators acting on 
the encoded qubit. (b) A manifestation of colored string as a 
bipartition. It stands for an error. 



and uiB disjoint regions, respectively. So, the cardinali- 
ties will be 

dA = 2^-*+2™s-2 ^ ^ 2^«+2™-*-2 . (40) 

Now we move on in order to calculate the topological 
entanglement entropy defined in Eq. (l38|) . For the parti- 
tions shown in Fig.(l4|) we have 

miB = TUiA = 2, 



miA = m2A = "^2_B = n^3A = W-3B = »7^4_B = 1, (41) 



SiA + = E2A -I- T.3A ■ (42) 

Finally the topological entanglement entropy reads as 
follows: 

Stopo = -4. . (43) 

Topological entanglement entropy depend only on 
the topology of the regions and no matter with their 
geometriesi^. This derivation for topological entan- 
glement entropy of color codes is consistent with the 
subleading term of scaling of the entropy. The topolog- 
ical entanglement entropy is related to the subleading 
term of the entanglement entropy, i.e. Stopo = —27 — . 
This leads to 7 = 2 which has also been derived Eq. ([32|) . 



V. ENTANGLEMENT PROPERTIES OF 
PLANAR COLOR CODES 

An important class of topological stabilizer codes in 
practice is planar codes which are topological codes that 
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can be embedded in a piece of planar surface. These pla- 
nar codes are very interesting for topological quantum 
memory and quantum computation o^'^^ . A planar color 
code will be obtained from a 2-colex without border when 
a couple of plaquettes removed from the code. For exam- 
ple when a plaquette, say green, is removed only green 
strings can have endpoint on the removed plaquette not 
the blue and red strings. The same scenario holds for 
other plaquettes. By this construction we will end with 
a 2-colex that has three borders each of one color. To 
each border only strings can have endpoint which have 
same color of that border. The most important class of 
such planar color code is triangular code that has been 
shown in Fig.®. The essential property of such code 
is determined via realizing a 3-string operator and its 
deformation. We denote them by T-^ and so that 
{T-^ ,T^} = since they cross each other once at strings 
with different colors. This latter anticommutation rela- 
tion can also be invoked by considering all qubits which 
makes the triangular color code very interesting for full 
implementation of the Clifford group without need for se- 
lective addressing^-. Instead of giving such fruitful prop- 
erties, we are interested in the entanglement properties 
of the triangular code. 

To have a concrete discussion, we consider a triangular 
lattice which contains hexagons like what has been shown 
in Fig.®. The total number of plaquettes and vertices 
are \P\ = |fc(fc + l) and \V\ = 3fc(fc + l)-M, respectively, 
where k in an integer number. The plaquette operators 
like those defined in Eq.([T|)have been used in order to 
define a stabilizer subspace for this code. A main point 
about the planar code is that there are no constraints on 
the action of plaquette operators, i.e they are indepen- 
dent. For example if we product all red and green pla- 
quette operators, the resulted operator will not be same 
with the product of red and blue plaquette operators, a 
property that was absent for color code on compact sur- 
face. However, the product of distinct plaquettes is in the 
stabilizer group. The coding space is spanned with the 
states which are fixed points of all plaquette operators. 
With the above identifications for the number of plaque- 
ttes and vertices, the dimension of the coding space will 
be A = ^ | p|^|_p = 2. This implies that triangular code 
encodes a single qubit. Two states are completely deter- 
mined by the elements of stabilizer group and 3-string op- 
erator. Note that here 3-string operator plays a role like 
nontrivial closed strings in the color code on the torus. 
For triangular code the 3-string is nontrivial, i.e. it com- 
mutes with all plaquette operators and is not product of 
some plaquettes. So the stabilized sates will be 

|0) = |G|-i/2^g|0)«'l^l , |1)=T^I0) . (44) 

A generic state can be in a superposition of the above 
states such as = X]i=o^*l*) where X]i=o ~ 
The reduced density matrix in Eq. (|lip gives what we 
need to evaluate the entanglement entropy. Again we 



make different bipartitions and calculate the entangle- 
ment entropy. 

Let the system be in the state |0). A single qubit is 
maximally entangled with rest of the system while two 
qubits are not entangled implying that the triangular 
code is also a multipartite entangled structure. As an 
another bipartition we consider a colored spin chain 
such as red chain in Fig.([5)D). This string contains 2k 
qubits and anticommutes with plaquettes which the 
string has an endpoint on them, i.e they share in odd 
number of qubits, thus its action on the coding space 
amounts to error or on the other hand to the appearance 
of excitations above the ground state if we adopt the 
coding space as ground state of a local Hamiltonian 
as in Eq. ()10p . However, such a string has maximal 
entanglement Sa — 2fc with other qubits of the lattice. 
This means that all configuration of the spin chain are 
allowed in the state |0) of the code, a feature that is 
similar to the entanglement of an open string in the 
topological color code on the torus (see Subsec(C)). 
Taking into account the qubits living on the 3-string net 
as subsystem A leads to occurrence of configurations of 
subsystem A with an only even number of spin flipped 
in the state |0). For this string net as shown in Fig. ((6^) 
the number of qubits of the string is 2fc + 1, and the 
entanglement entropy then will become Sa — 2fc. By 
inspection we see that the plaquette operators flip only 
even number of qubits of the string net. Now let the 
system be in a generic state such as All possible 

configurations of string net are allowed since all of them 
can be realized by applying some 3-string operators T-^ . 
There are 2^*"' configurations with even number of spins 
flipped and 2^*^ configurations with odd number of spins 
flipped. The states |Q) and |l)include even and odd 
configurations, respectively. Therefore the entanglement 
entropy for a generic state will be: Sa = 2k + H{a) 
where a — |ao|^. 

What about the topological entanglement entropy of 
planar color codes? In order to answer this question 
we obtain a set of bipartitions such as in Fig.® in a 
large triangular code. Let and stand for the 
number of plaquette operators acting only on A and 
-B, respectively, and Sab stands for the number of 
plaquette operators acting simultaneously on A and 
B. All plaquette operators are independent. To get 
things simpler, let consider a simple case where the 
subsystem A is a single convex connected region. The 
cardinality of the subgroup Ga is 2^^. However, there 
are two non-trivial closed string acting on subsystem 
B which is impossible to provide them by product of 
some plaquettes of -B. In fact these two independent 
strings are resulted from the product of some plaquettes, 
say red and green, which are free to act on A. So the 
cardinality of the subgroup B reads 2^^+^. Thus for the 
entanglement entropy we obtain Sa = '^ab — 2. 
In order to calculate the topological entanglement 
entropy, first the entanglement entropy is calculated for 
different bipartitions and then the expression in Ea. ((55|) 
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gives the topological entropy. Because of the boundaries 
in the planar codes, for a disconnected region like (1) or 
(4) in Fig.Q the constraint we imposed in Eq. (|40|) is no 
longer true. For the bipartition (1) in Fig.Q the cardi- 
nalities of subgroups Ga and Gb will be c?^ — 2^^^'^'^ 
and ds = 2^^^+^, respectively. However, for the bipar- 
titions (4) in Fig.lg]) the cardinalities wiU be = Z^''-^ 
and ds = 2^-*^+*, respectively. Thus the topological 
entanglement entropy becomes Sf^p^ — —4. Although 
the closed string structures of bipartitions in the tri- 
angular code due to the boundary effects differ from 
the color code on the torus, the topological entangle- 
ment entropy is same for both structures that is related 
to the fact that both structures have the same symmetry. 



VI. CONCLUSIONS 

In this paper we calculated the entanglement prop- 
erties of topological color codes (TCC) defined on a 
two-dimensional lattice. The entanglement entropy 
was measured by the von Neumann entropy. We 
considered two structures of TCC either defined on the 
compact or planar surface. The coding space of TCC 
is spanned by a set of states in which are the fixed 
points of a set of commuting Pauli operators, i.e they 
are common eigenvectors of all elements of stabilizer 
group with eigenvalue +1. In fact this set stabilizes 
the coding space. The stabilizer group is generated 
from plaquette operators. Product of different plaquette 
operators produces colored strings which in fact are the 
closed boundary of several plaquettes. However, There 
exist some independent non-trivial strings winding the 
handles of manifold where the model defined on. This 
non-trivial strings commute with all plaquette operators 
but are not actually in the stabilizer group. These 
nontrivial strings make remarkable properties of the 
color code more pronounced. 

For a manifold with genus g, the coding space spanned 
by 4^8 states that can be adopted in which be ground 
state of a local Hamiltonian. The degeneracy of ground 
state subspace depends on the genus of the manifold 
which is a feature of topological order. Different states 
of ground state subspace can be constructed by means of 
spin flipped elements of stabilizer group and nontrivial 
closed loops. In order to calculate the entanglement 
properties of the TCC, the reduced density matrix of a 
subsystem is obtained by integrating out the remanding 
degrees of freedom. We did this by using the group 
properties of the stabilizer group. 

For both structures of color code on compact and planar 
surface while a single qubit is maximally entangled with 
others, two qubits are no longer entangled. This finding 
manifests the color code is a genuinely multipartite 
entangled structure. We also considered other bipar- 
titions such as spin chains, spin ladders with various 
colors and homologies. For all bipartitions we found the 



entanglement entropy depends on the degrees of freedom 
living on the boundary of two subsystems. However, in 
the entanglement entropy of a convex region, there is 
a subleading term which is ascribed to the topological 
properties of region not the geometry, a feature that has 
also been arisen in a massive topological theoryiS. The 
fact that the entanglement entropy for a region of lattice 
scales with its boundary is a feature of area law which 
is also of great interest in other branch of physics such 
as black holes^^. 

We exploited the scaling of the entropy to construct a 
set of bipartitions in order to calculate the topological 
entanglement entropy of color codes. We find that it 
is twice than the topological entropy in the topological 
toric code. The non-vanishing value for topological 
entanglement entropy and dependency of degeneracy 
on the genus of the surface demonstrate remarkable 
features on the fact that topological color code maybe 
fabric to show topological order. For the toric code 
model the total quantum dimension of excitation is 4 
which stands for different superselection sectors of its 
anyonic excitations. The total quantum dimension of 
color code is 16 inspiring that the topological color 
codes may carry richer anyonic quasiparticles. They 
may be colored and their braiding will have nontrivial 
effect on the wavefunction of the system. Indeed 
braiding of a colored quasiparticle around the another 
one with different color will give rise to global phase 
for the wavefunction, i.e they are abelian excitations. 
Triangular color code with a remarkable application for 
entanglement distillation is also a highly entangled code. 
Although the colored strings employed in the 2-colex in 
the compact surface are not relevant in the triangular 
code, the notion of string-net gives some essential 
properties to the triangular code. In the state of the 
triangular code all configurations of a colored string, say 
a red one in Fig.([5}3), are allowed, i.e is has maximal 
entanglement with the rest of the system. However, for 
a string-net only the configurations with even number 
of spin flipped are allowed. Entanglement entropy scales 
with boundary of the region and incudes a topological 
term. This topological terms yields a non-vanishing 
topological entanglement entropy which is same with 
what we obtained for the color codes on the torus. 
Similar entanglement properties of color codes defined 
on the compact and planar surface are not a matter of 
chance because both structures carry same symmetry. 



I would like to thank A. Langari for fruitful discussions 
and comments. I would also like to acknowledge M. 
A. Martin-Delgado for his useful comments. This work 
was supported in part by the Center of Excellence in 
Complex Systems and Condensed Matter (www.cscm.ir) . 



12 



References 



* 'kargarian@physics.sharif.edu 

^ S. Oh and J. Kim, Phys. Rev. B 71, 144523 (2005). 

2 V. Vedral, New J. Phys 6, 102 (2004). 

^ X. G. Wen, Phys. Lett. A 300, 175 (2002). 

* L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 
Phys. 80, 517 (2008). 

^ N. Goldenfeld, lectures on phase transitions and the renor- 

malization group, Westview Press, 1992. 
" X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990); 

X.G.Wen, Phys. Rev. Lett. 90, 016803 (2003). 
^ A. Y. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 

* D. Gottesman, Phys. Rev. A. 54, 1862 (1996). 

^ E. Dennis, A. Kitaev, A. Landahl and J. Preskill, J. Math. 
Phys. 43, 4452 (2002). 

A. Hamma, R. lonicioiu, and P. Zanardi, Phys. Rev. A. 
71,022315 (2005); A. Hamma, R. lonicioiu, and P. Zanardi, 

Phys. Lett. A 337, 22 (2005)^ 

" S. Ryu and T. Takayanagi,'arXiv:hep-th/0605073 

" A. Kitaev and J. PreskiU, Phys. Rev. Lett. 96, 110404 



(2006). 

M. Levin and X. G. Wen, Phys. Rev. Lett. 96, 110405 
(2006). 

H. Bombin, M. A. Martin-Delgado, Phys. Rev. lett, 97, 
180501 (2006). 

H. Bombin, M. A. Martin-Delgado, Phys. Rev. A, 76, 
012305 (2007). 

H. Bombin, M. A. Martin-Delgado, Phys. Rev. B, 75, 
075103 (2007). 

^"^ R. Raussendorf, S. Bravyi, J. Harrington, Phys. Rev. A, 
71, 062313 (2005). 

W. K. Wootters, Phys. Rev. lett, 80, 2245 (1998). 
M. B. Plenio, J. Eisert, J. Dreifiig, and M. Cramer, 
Phys.Rev.Lett. 94, 060503 (2005). 
^° R. Bousso, Rev. Mod. Phys. 74, 825 (2002). 

S. Ryu and T. Takayanagi, Phys. Rev. Lett. 96, 181602 
(2006). 



